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Abstract 

We present the results of an analysis of data collected by IceCube/DeepCore in 2010-2011 resulting in the first 
significant detection of neutrino oscillations in a high-energy neutrino telescope. A low-energy muon neutrino sample 
(20-100 GeV) containing the oscillation signal was extracted from data collected by DeepCore. A high-energy muon 
neutrino sample (100 GeV -10 TeV) was extracted from IceCube data in order to constrain the systematic uncertainties. 
The non-oscillation hypothesis was rejected with more than 5cr. We fitted the oscillation parameters Am^ and sin 2 2023 
to these data samples. In a 2-flavor formalism we find Am 2 3 = (2.5 + 0.6) • 10~ 3 eV 2 and sin 2 2023 > 0.92 while 
maximum mixing is favored. These results are in good agreement with the world average values. 
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1. Introduction 

Neutrino oscillation experiments have established 
that neutrino flavor and mass eigen states do mixJT|. 
So far, solar and long-baseline reactor neutrino ex- 
periments have measured the mass-mixing parameters 
(5m 2 , 812) in the v e — > v e channel (electron neutrino dis- 
appearance), while atmospheric and long-baseline ac- 
celerator experiments have measured (Am 2 , #23) in the 
v pt - > Vfi channel (muon neutrino disappearance^] 

We use data collected from May 2010 to May 2011 
by the IceCube neutrino telescope with its low-energy 
sub-detector DeepCore |2| to measure the atmospheric 
neutrino oscillation parameters. The IceCube Neutrino 
Observatory is a cubic -kilometer neutrino detector in- 
stalled in the ice at the geographic South Pole Q. It 
is based on the optical detection of secondary parti- 
cles produced by neutrinos interacting in the ice or the 
bedrock below. Such charged particles emit Cherenkov 
light, which is detected by IceCube's optical sensors. 
These sensors are attached to 86 strings, which hold 60 



1 We here adopt a convention where dm 2 = m 2 - m 2 and 
Am 2 = m 2 - (m^ + m 2 )/2. 



sensors each. This corresponds to a standard vertical 
spacing of 17 m between sensors and a horizontal dis- 
tance of 125 m between the strings. The DeepCore sub- 
detector consists of eight additional strings deployed 
in the center of IceCube and the surrounding IceCube 
strings. On the dedicated DeepCore strings, the sensors 
are concentrated in the cleanest deep ice, resulting in a 
denser 7 m vertical spacing of sensors there. During the 
data taking period of this analysis, 79 detector strings 
were operational (IceCube-79) including six dedicated 
DeepCore strings. 

We used a 2-flavor formalism to describe neutrino os- 
cillations, neglecting 3-flavor effects and matter effects. 
In this formalism, the muon neutrino survival probabil- 
ity is given by 

P(y M -+v M ) = l- sin 2 26»23 sin 2 (1.27Am 2 L/£') (1) 

with L as the length of propagation of the neutrino in 
km and E as the neutrino energy in GeV. 

2. Data sample 

We extracted two samples of upwards going neutrino 
events from data collected by IceCube-79, one at rela- 
tively high energies using data from the entire IceCube 
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Figure 1 : Distribution of the neutrino energy of atmospheric neutri- 
nos in the low-energy (DeepCore) and in the high-energy (IceCube) 
sample. 



Figure 2: Zenith resolution (mean absolute difference between true 
and reconstructed track) for the low-energy sample as a function of 
true zenith. 



detector and one at lower energies selected in the Deep- 
Core volume, rejecting backgrounds by an active veto 
|4|. Neutrino oscillations are expected to affect only the 
low-energy sample. The high-energy sample provided 
large statistics outside the signal region and served to 
constrain systematic uncertainties. 

The directions of the muon tracks in the high-energy 
sample were reconstructed with the standard maximum 
likelihood muon track in IceCube [5|. For low-energy 
events, the same method was applied as an initial step. 
As the hypothesis of a throughgoing track is not correct 
for these, the finiteness of the tracks was considered in 
a second step. The track length and the start/stop points 
of the track were determined as well as the likelihood 
whether the track is starting and/or stopping in the de- 
tector [4]. Quality cuts like the number of unscattered 
photons and the track likelihood allowed for the rejec- 
tion of misreconstructed downwards going muon back- 
ground. 

In Fig. [T] the neutrino energy distribution of the low- 
energy and the high-energy sample are shown. The 
resolution of the reconstructed zenith angle is essen- 
tial because the propagation length is proportional to 
cos(zenith). A variation of the zenith thus represents a 
variation of L/E. As displayed in Fig. [2] a resolution of 
8° is achieved for the low-energy sample, independent 
from the zenith. 

A modified version of the atmospheric neutrino flux 
model derived by Honda et al. [6| was used because re- 
cent measurements of the spectrum of charged cosmic 
rays in the energy range of 200 GeV to 100 TeV indi- 
cate a cosmic ray spectrum flatter than that assumed by 
Honda et al., see e.g. (7j. We adjusted the resulting neu- 
trino spectrum by hardening its spectral index by 0.05 to 
reflect these new measurements. 



3. Systematic uncertainties 

A covariance matrix in a x 2 fit was use d to consider 
systematic uncertainties in the data analysis. In order to 
obtain the most likely value of the individual sources of 
systematic errors, the pulls as defined in [8] were used. 
The following sources of systematic uncertainties were 
considered explicitly and propagated by Monte Carlo 
(MC) simulation to the final selection level: 

• the absolute sensitivity of the IceCube sensors 
(±10%) and the relative efficiency of the more effi- 
cient DeepCore sensors (1.35 ± 0.03) 

• the optical parameters (scattering, absorption) of 
the ice as detector medium: the uncertainty is es- 
timated by the difference of the optical parameters 
obtained by the extraction methods [9| and IfTUl 

• the absolute normalization of the cosmic ray flux 
(+25%) and its spectral index (±0.05) 

• the uncertainty of the neutrino production rate in 
the atmosphere: the difference of calculations by 
[6 1 and [ 1 1 1 were used for and for v e . 

4. Results 

During May 2010 to May 2011, we collected 318.9 
days of high quality data, excluding periods of cali- 
bration runs, partial detector configurations and detec- 
tor downtime. The low energy sample contained 719 
events, while the high energy sample contained 39, 638 
events after final cuts. In a first step, we evaluated the 
X 2 for the data collected by IceCube for two different 
physics hypotheses: the standard oscillation scenario 
represented by the world average best fit parameters and 
the non-oscillation case. With A^ 2 = 30, we rejected 
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Figure 3: 68% and 90% CL contour line of the result of the IceCube- 
79 oscillation analysis in comparison with the results of Antares| 12], 
Minos and SuperKamiokande l 131 . 

the non-oscillation hypothesis with a p-value of 10~ 8 , 
corresponding to 5.6<x. The significance was evaluated 
by a toy MC considering deviations from a x 2 distri- 
bution (both hypotheses do not correspond to the mini- 
mum^ 2 ). 

In a second step, the x 2 was evaluated as a function 
of the oscillation parameters. The best fit is given by 
Am?, = 2.5 ■ 10~ 3 eV and sin 2 26*23 = 1, with an absolute 
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Figure 4: Data and MC expectation at best fit (physics parameters and 
pulls) in the low energy sample and for the high energy sample. 



X 2 = 11.3 and 18 degrees of freedom (20 bins, 2 fitted 
parameters). The value of the absolute x 2 corresponds 
to a goodness-of-fit p-value of 0.88. This indicates a 
good agreement of data to MC within the assumed un- 
certainties. All pulls on the systematic uncertainties are 
within the l<x uncertainty range. The data as a func- 
tion of zenith together with the statistical lcr range of 
the MC expectation corrected for the pulls is shown in 
Fig. |4] The result is in good agreement with other ex- 
periments, which measured the atmospheric oscillations 
with a high resolution at lower energies! 1 3 1 . 

The two dimensional confidence regions of the os- 
cillation parameters in this measurement were deter- 
mined from the Ay 2 around the best fit with two de- 
grees of freedom (68% CL: A^ 2 = 2.30 and 90% CL: 
A^ 2 = 4.61), see Fig. [3] The confidence regions of the 
individual parameters were determined by marginaliza- 
tion analogous to a profile likelihood method. We obtain 
68% CL intervals of Am 2 3 = (2.5 + 0.6) • 10~ 3 eV 2 and 
sin 2 26*23 > 0.92 using Ax 2 with one degree of freedom. 

The analysis of IceCube data presented here provided 
the first significant detection of atmospheric neutrino 
oscillations with a high-energy neutrino telescope. In 
future, a significant improvement of the resolution of 
IceCube on the atmospheric neutrino oscillations is ex- 
pected by the inclusion of the reconstructed neutrino en- 
ergy in the analysis, by the use of new reconstruction 
methods which are more efficient at lower energies and 
by the inclusion of two additional dedicated DeepCore 
strings which started data taking in May 201 1. 
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